Disorder from disorder in a strongly frustrated transverse field Ising chain 
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We study a one-dimensional chain of corner-sharing triangles with antiferromagnetic Ising inter- 
actions along its bonds. Classically, this system is highly frustrated with an extensive entropy at 
T — and exponentially decaying spin correlations. We show that the introduction of a quantum 
dynamics via a transverse magnetic field removes the entropy and opens a gap, but leaves the ground 
state disordered at all values of the transverse field, thereby providing an analog of the "disorder 
by disorder" scenario first proposed by Anderson and Fazekas in their search for resonating valence 
bond states. Our conclusion relies on exact diagonalization calculations as well as on the analysis 
of a 14th order series expansion about the large transverse field limit. This test suggests that the 
series method could be used to search for other instances of quantum disordered states in frustrated 
transverse field magnets in higher dimensions. 

PACS numbers: 05.50.+q, 75.10.-b, 75.10.Jm, 75.30.Kz 



I. INTRODUCTION 

The study of quantum versions of classically frustrated 
magnets has been a subject of interest at least since 
the work of Anderson and Fazekas on the possibility 
of a quantum disordered state for the triangular lattice 
Heisenberg antiferromagnetHI Following their work, and 
especially after the early suggestion of Anderson that the 
cuprate superconductors derive their special properties 
from the proximity of a spin liquid stateJ3 there has been 
a considerable amount of work on quantum Heisenberg 
models on various frustrated lattices such as the kagome.0 

One can think of quantum Heisenberg models as (clas- 
sical) Ising models perturbed by a transverse (XY) 
exchange — indeed, this was the strategy followed by An- 
derson and Fazekas in their analysis of the triangular 
lattice problem. Phrased in this fashion, the problem be- 
comes one of the effects of introducing a quantum dynam- 
ics into a highly degenerate ground state manifold of the 
Ising system — a procedure of evident interest on account 
of the singular effects of any perturbation. The simplest 
instance of this more general problem is the introduc- 
tion of a transverse magnetic field, which has been used 
in several contexts previously to generate the simplest 
quantum statistical mechanics; for a review see Ref. |] 
which also reviews some work on some one-dimensional 
systems with competing interactions. 

A number of geometrically frustrated quantum Ising 
systems have been studied recently by Moessner, Chan- 
dra and one of the present authors who have reported in- 
stances of "order by disorder"0 in which a non-trivial or- 
dering pattern is selected by the quantum fluctuations as 
well as one of "disorder by disorder" (on the kagome lat- 
tice) in which the Anderson-Fazekas scenario of a disor- 
dered quantum state constructed out of a disordered clas- 
sical manifold is realized.19 In this paper, we report stud- 
ies of a one-dimensional, geometrically frustrated chain 
which also exhibits disorder by disorder. This is interest- 
ing in itself; additionally, through a technical advance, we 



were able to provide strong evidence for the quantum dis- 
ordered state by means of a strong coupling perturbation 
expansion, which we anticipate will be a useful technique 
in studying higher-dimensional frustrated quantum Ising 
systems. 



II. THE MODEL AND ITS POSSIBLE 
ORDERING 



We consider the Hamiltonian 

W) 1 



(1) 



where J > is the antiferromagnetic Ising exchange, the 
sum in the first term runs over the bonds of the triangular 
chain in Figure |l, T is the strength of the transverse field 
and the a a are the Pauli operators for S = 1/2. (As the 
true spin operators are S a — cr a /2, when U = 1, our J 
is really one quarter the exchange and our T is half the 
physical transverse field.) 
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FIG. 1. The triangular Ising spin chain. 
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FIG. 2. A typical ground state at T — 0. The dashed lines 
represent frustrated or "bad" bonds, of which there is exactly 
one per triangle. 



The ground states of the system without the trans- 
verse field can be obtained by minimizing the energy of 
each triangle separately and this in turn requires that we 
choose the one bond in each that will remain unsatisfied 
(the "bad" bond). The choice of the bad bonds is totally 
independent from triangle to triangle (the system is a bad 
bond paramagnet); a typical state is exhibited in Fig ||. 
With free boundary conditions, the number of such states 
is thus 2x3* where N is the number of triangles and 
the extra factor of 2 accounts for the two Ising reversed 
states that give rise to the same bond configuration.!!! 
This extensive entropy is accompanied by short-ranged 
correlations and we find that the spin-spin correlation 
function averaged over the ground state manifold takes 
the form 

(S z at Sf jt ) = (S^Sfa) = (-l/3f~ a+1 (2) 
(S' ab S' pb ) = (S* at S% b ) = (-l/'Sf- a (3) 

where at and ah denote the top and bottom spins in 
unit cell a. Due to the lack of inversion symmetry in 
our choice of unit cell, these forms hold for [3 > a; they 
can be extended to cover (3 < a by inspecting the lat- 
tice and noting its symmetries. Evidently, the system is 
disordered at all temperatures. 

We turn now to the quantum problem posed by the 
inclusion of the transverse field. 



A. Analysis at V < J 

At T — 0, an infinitesimally small T is a singular per- 
turbation which will lift the macroscopic degeneracy of 
the classical problem. Following Ref. ^, we will attempt 
to identify the physics of this regime variationally. To 
this end we identify the configuration (unique, up to 
global symmetry operations) that maximizes the num- 
ber of "flippable" spins, which is shown in Figure || A 
flippable spin in a given configuration is one which can 
be reversed without violating the ground state constraint. 
As a flippable spin can be polarized along x at no cost in 
exchange, we expect a state that maximizes their number 
to be especially favored by the transverse field. As not all 
flippable spins are independently flippable, we have two 
options. We may construct the staggered state in which 
we polarize the maximal set of independently flippable 
spins along the top row, 

|s) = (®to P |a:-» | UUU • • ^bottom ; (4) 



here \x— ) is the state with a spin pointing down along the 
x axis. Alternatively we may construct the uniform state 
in which we polarize all the flippable spins but correct for 
the conflicts by projecting onto the ground state manifold 
by the action of the projector P which eliminates those 
configurations that do not have exactly one bad bond per 
triangle: 

|«> =^{(®to P |z-» (5) 

® (| T)<8 ® | t)<8 \X-) ® | t) • • Obottom} ■ 

The uniform state inherits all the symmetries of the maxi- 
mally flippable configuration but the staggered state does 
not. For the wave functions as written, we can readily 
evaluate the energies and we find that the staggered state 
has the same energy — (r + J)/2 per site as the uniform 
state. We can also write down the magnetization (up to 
translations and Ising reversal) in the staggered state: 

and the uniform state 

m{a) = {~f+(-lr bottom. ^ 

Needless to say, this analysis being purely variational is 
only indicative of what kind of "order from disorder" the 
system might exhibit at T <C J- At issue is whether 
the connectivity of the set of configurations is too high 
(e.g. as reflected in the disorder in the classical prob- 
lem) to permit a localization of the ground state wave- 
function near the maximally flippable configuration. As 
advertised we will, instead, find that the system remains 
disordered. 

AAAAAAAA 

FIG. 3. The maximally flippable configuration. All sites 
along the top row and half of those along the bottom row are 
flippable. 



B. Analysis at T > J 

One of the attractive features of transverse field Ising 
models, emphasized in Ref. [| is that they exhibit a 
gapped, paramagnetic phase at T ^> J. This enables 
a natural perturbative expansion for the problem, in 
which we perturb in the Ising exchange about the purely 
transverse field paramagnet. A second attractive feature- 
exhibited by the well known Trotter-Suzuki procedures 
is that the quantum partition function (with the trans- 
verse field present) has the form of the classical partition 
function of copies of the chain (now without the trans- 
verse field) coupled ferromagnetically in the imaginary 



2 



time direction which can then be analyzed from its "high 
temperature" phase via the standard Landau-Ginzburg- 
Wilson (LGW) procedures It is not difficult to show that 
these are attempts to search for the same instability. We 
will pursue the perturbative expansion a bit later; here, 
we will attempt to identify the leading instability from 
the paramagnetic phase by the LGW step of diagonaliz- 
ing the bond (adjacency) matrix. (The equivalence to the 
first step in the perturbation expansion will be evident 
when we come to it.) 

As the ferromagnetic interaction in the imaginary 
time direction is trivially incorporated by including a 
quadratic dispersion about k T = 0, we focus on the prob- 
lem in the chains. To diagonalize the quadratic form, 



mean-field estimate for the location of the transition. To 
this end we consider the possibility of ordering into the 
staggered state. In this case only the spins on the bottom 
row sec an effective field and the problem (in mean field 
theory) is identical to that of the purely one-dimensional 
transverse field Ising chain. This yields a critical coupling 



x c = (J/r) c = 1/2 

and a staggered magnetization 



W = zl S iSj 



(8) 



we introduce Fourier variables with respect to the unit 
cell location, 



Sat — ^ e lka Skt 



Sab — ^ 



Aka n 



(9) 
(10) 



and find that the Hamiltonian takes the form 



w = ^2[ S H s kb ] 



k>0 



(1 



-ik\ 



(1 



ik\ 



2 cos k 



Skt 

Skb 



(11) 



Further diagonalization yields the eigenvalues 



A± = (e - 1) ± Ve 2 + 1 



(121 



where we have defined e = 1 + cos k. Evidently, A_ is 
minimized when one has e = — > k = tt. The corre- 
sponding eigenvector is 



(13) 



Hence, the spatial dependence of the lowest eigenmode, 
which is the candidate for ordering, is 



(-l) c 



(14) 



Thus, as r/J is lowered — which corresponds to lower- 
ing the temperature in the classical representation — the 
LGW analysiSpSuggests a transition in the d = 2 Ising uni- 
versality classli3 to a state with the top row disordered 
and bottom row ordered antiferromagnctically, which is 
the staggered state that we had constructed earlier. 



C. Mean Field Theory 

As we are interested in the suppression of ordering by 
fluctuations due to the frustration, it is useful to have a 



y/l ~ (2z)" 2 



(15) 



(16) 



in the ordered phase. The magnetization of the top row is 
always zero. Better estimates can be obtained by treating 
the top row in mean field theory .and using the known 
exact results for the bottom row.Eil In this fashion we 
obtain x c = 1 and m ~ (x — ic) 1 / 8 . 



III. EXACT DIAGONALIZATION 
A. Modified Lanczos method 

We now turn to our numerical studies of the model, 
beginning with exact diagonalization of the Hamiltonian 
for finite lattices with periodic boundary conditions. Ex- 
act diagonalization unfortunately imposes a demand on 
computer memory which rises exponentially with the size 
of the system considered. Consequently, the systems dis- 
cussed in this report are limited to at most 8 unit cells 
(16 spins). _ 

In the modified Lanczos method that we applyjLj one 
begins with a state |0) expected to have nonzero overlap 
with the ground state of 7i. For the present problem we 
used the trivial (fully x-polarized) J = ground state, 
which is a q = state. One generates a state orthogonal 
to |0) by acting upon |0) with H and using Gram-Schmidt 
orthogonalization to generate the state 



|1> = W|0>- 



(0|0) 



(17) 



One then diagonalizes TL in the resulting two-dimensional 
subspace. The lowest eigenvalue is the improved estimate 
for the ground state energy, and the corresponding eigen- 
vector is the improved estimate for the ground state wave 
function. This process is iterated to convergence. 

A slight modification of the procedure used in obtain- 
ing the ground state energy and wave function can be 
employed to find the energy gap and first excited state 
wave function. One begins with an initial guess with 
nonzero overlap with the lowest excited state wave func- 
tion and orthogonal to the ground state, and then pro- 
ceeds as described above for the ground state. The initial 
guess for the excited state was chosen to be one of the 
excited states in the J = system, with a single flipped 
spin at a particular location. Although the state is not 
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one of definite momentum, many iterations of the modi- 
fied Lanczos procedure converge to a state of definite mo- 
mentum in which k = tt, as expected from the analysis in 
Section 0. Although successive excited state estimates 
should in principle be orthogonal to the ground state, 
rounding errors cause an admixture with the ground state 
to occur at each iteration. To prevent a convergence to- 
ward the ground state, Gram-Schmidt orthogonalization 
is used at each step to remove any component of the 
ground state which may have entered. Although one does 
not have the exact ground state at a particular iteration, 
convergence can still be achieved by performing Gram- 
Schmidt orthogonalization with the most recent estimate 
for the ground state wave function. 



B. Excitation gap 

Figure || displays the gap to the first excited state com- 
puted for a variety of system sizes (N denoting the num- 
ber of unit cells) and values of x = J/T. For convenience, 
the domain [0, oo] of x is compressed into the interval 
[0, f] via the transformation y = x/(x + 1). For most val- 
ues of y, convergence of the gap with increasing system 
size is quite rapid, in itself an indication of the presence 
of a short correlation length. In fact, for y < 2/3, the 
gaps for systems with 6 and 8 unit cells are virtually 
indistinguishable. In that regime, it is clear that the sys- 
tem remains disordered (in the same phase as the pure 
transverse- field model y — 0). For larger y the situation 
is less clear-cut, so we next consider a more sensitive test 
for the existence of a continuous phase transition. (We 
note that our largest system size yields a ground state 
energy of — \ [J + (l.ll)r] per spin at T <C J, 11% lower 
than the energy for the variational states considered ear- 
lier.) 
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FIG. 4. Gaps obtained via exact diagonalization for 
N = 2,4,6,8 



C. Cumulants 



Following Binder ,0 we examine the cumulants 

_ d*l 4 ) 



(18) 



In equation |i~8|, ^ is the order parameter of interest. Since 
any ordering is expected to be accompanied by antiferro- 
magnetic ordering of the bottom row, a sensible choice of 
\1/ seems to be the staggered magnetic moment operator, 



M, 



Mr 



(19) 



on the bottom row. 

It is expected that for systems of size L, in the vicinity 
of a critical point, r obeys the scaling form 



= f((x-x c )L 1 /») 



(20) 



If a critical point exists, the cumulants computed for var- 
ious system sizes should intersect at the same point for 

Cumulants for systems containing 4, 6, and 8 unit cells 
are displayed in Figure ^|. There is apparently no in- 
tersection of the cumulants for different system sizes for 
x < 20 (and, by extrapolation, for even larger values of 
x). The behavior of the cumulants is consistent with the 
absence of a critical point and therefore suggests that the 
system is disordered for all values of J/T. We note that 
in a disordered phase with a Gaussian distributed order 
parameter, the cumulant will approach the value 3. At 
x = the cumulant can be shown analytically to have 
the value 3 — 4/iV for N unit cells; evidently the finite size 
effects are greater at large x, where the ground state con- 
straint produces correlations, but even there the steady 
growth with system size is consistent with a thermody- 
namic limit value of 3. 



y = x/(x+1);x = J/T 
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FIG. 5. Plots of the cumulant (M S 4 )/|(M S 2 )| 2 for triangular 
chains containing 4, 6, and 8 unit cells. 



IV. PERTURBATION EXPANSION 
GENERATION AND ANALYSIS 

A. Cluster expansions for the elementary excitation 
spectra 

A perturbation expansion for the energy gap (and the 
full elementary excitation spectrum, incidentally) about 
the J = limit was constructed to 14th order in J/Y 
using a cluster expansion technique. The ideas under- 
lying this calculation are set out briefly in Ref. [l4|, and 
described in more detail in Ref. [D| However, this partic- 
ular calculation-is unusual (but not unique, see Muller- 
Hartmann et a/Jl3) in one interesting respect: it is carried 
out for a system with more than one branch of elemen- 
tary excitations which are degenerate in the unperturbed 
limit, but are not degenerate by symmetry and hence are 
not degenerate for J > 0. That such calculations can be 
done by cluster expansion methods was mentioned briefly 
in Ref. [l5|. Here we supply some further details which are 
important in making the connection between the the im- 
mediate product of the cluster expansion, which is a ma- 
trix containing power expansions that describe the mo- 
tion of elementary excitations in real space, and the de- 
sired end product, the power series expansions for the 
excitation spectra in wave vector space. 

For the purposes of this section we assign the sites at 
the base of the triangles the coordinates (n, 0), and those 
at the tips of the triangles (n+ 1, 1), with n running over 
the integers. The unperturbed Hamiltonian is the trans- 



verse field J2nl a 



(n,0) 



V+i,i) 



and the perturbation is 



the Ising exchange £ n [ ff (n,o) tr (n+i,o) + *(»,(»*(«+$, i) 
a ( n +i o) a (n+i i)]' ^ ue P er turbation expansion is a power 
series in x = J/Y. 

For the unperturbed system, the ground state consists 
of all spins "down" (pointing along —x) and the elemen- 
tary excitations are single spins which have been flipped 
"up" . Note that there are two degenerate branches of 
such excitations, flat bands at ui = 2. (One could con- 
sider an unperturbed Hamiltonian of the more general 



form £ 



'(n,0) 



(n,l) 



which has the same symme- 
try as the triangular chain but lacks the degeneracy. For 
C 7^ 1, then, perturbations expansions for each of the 
branches of excitations could be constructed separately; 
however, there would be energy denominators (1 — C) in 
abundance that would lead to poorly converging series 
for £ near 1.) 

In our calculations, with £ = 1, we treat all of the 
elementary excitations on equal footing. We keep track 
of three distinct types of matrix elements in the effec- 
tive Hamiltonian, those that connect two base sites (type 
A), those that connect two tip sites (type B), and those 
that connect base to tip sites (type C); let us call t* 
the matrix element (which is a series expansion in A) in 
the effective Hamiltonian of type X that couples sites at 
horizontal distance d. The effective Hamiltonian has the 
following structure near the diagonal 



t c 



1/2 



L 3/2 



Z 5/2 



Z l/2 


l l 


b 3/2 
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l l/2 
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l l 


t C 

L 3/2 
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l l 


hl2 


]L 2 
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t C 
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t A 

H 


t C 


t A 

'o 


t C 


l 2 


L 3/2 


\g 
l l 


l l/2 


l 



V 



(21) 



This matrix acts on a column- vector describing the spin- 
flip amplitude at each site, listed in order of its x coordi- 
nate. This matrix is readily diagonalized by plane waves 
with a two-site basis. One obtains the two branches of 
the excitation spectrum in the form 



Fa + Fb±\{F a -F b 



F 2 c 



with 



F A(q) = ^2 t A cos nq , 

n>0 

likewise for Fs(q), and 

Fc(q) = f (2„+i)/2 cos(2n + l)q/2 

n>0 



(22) 



(23) 



(24) 



The relationship between the ts and ts is as follows. Take 
any two adjacent rows of the effective Hamiltonian, and 
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count how many tg there are for a given X and d: that 
number is the ratio t*/tf. Thus t£ = tfi, ?f /2 = 4tf /2 , 

and so forth. It is the is that come directly out of our 
computer programs and so it is convenient to express the 
functions that appear in Eq. ( p^ ) in those terms. 

We should note that our 14th order calculation of the 
excitation spectrum involved weight calculations for 1355 
graphs of up to 15 spins. We did not attempt to clas- 
sify the graphs topologically. The implementation of the 
weight calculations traded off considerable efficiency for 
safety (guarding against a variety of possible coding er- 
rors), and the calculations were carried out on a modest 
machine (333 MHz Pentium II with 256 megabytes of 
RAM), so it should be feasible to evaluate several more 
terms of the perturbation expansion if desired. 

B. Analysis of the energy gap series 

One is faced immediately with an interesting choice 
in analyzing the series for the gap, e_. One can either 
extrapolate the series for Fa, Fb and Fc individually 
and insert the results into the formula (|2|), or evaluate 
that formula for the series and extrapolate e_ directly. 

However, the lowest excited state is always found at 
q = tt, and in this case the need to choose vanishes be- 
cause Fc(ir) = and thus e_(?r) = min(i 7 ' J 4(7r), Fb(tt)) = 
F A (n). In Table J we display F A (n) and F B (n). The 
complete set of are available from the authors on re- 
quest. 



TABLE I. 


Coefficients of x n in the series Fa 


7r) and F b (tv). 


n 


Fa(tv) 


F b (tv) 





2.000000000 


2.000000000 


1 


-2.000000000 


2.000000000 


2 


1.000000000 


1.000000000 


3 


0.000000000 


-0.750000000 


4 


-0.250000000 


0.000000000 


5 


0.093750000 


-0.328125000 


6 


-0.144531250 


0.785156250 


7 


0.396484375 


-0.621093750 


8 


-0.544921875 


0.763671875 


9 


0.632377624 


-1.323387146 


10 


-0.908432960 


1.746130307 


11 


1.405481246 


-2.362823632 


12 


-2.073822043 


3.583402837 


13 


3.127851625 


-5.491168599 


14 


-4.983237137 


8.550544610 



To analyze the gap series, we applied the transforma- 
tion y = x/(x + 1), and constructed Pade approximants 
to the transformed series. We use the notation P^j, to de- 
note an approximant with a polynomial of order M in the 
numerator and order M' in the denominator. Figure ^ 
displays several approximants which use all the terms of 
the series (and hence M + M' = 14). 




-0.5 1 1 1 1 1 1 1 1 1 1 
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y = x/(x+1);x = J/T 

FIG. 6. Pade approximants and exact diagonalization 
results (TV = 8) for the gap plotted as a function of 
y = x/(x + l). 



A salient feature of Figure || is the consistency of the 
Pade approximants over a large fraction of the domain. 
The approximants displayed in the plot are consistent to 
within 1% for y as large as 0.75. Note that this corre- 
sponds to a value of a; = 3, well in excess of the mean field 
estimate of x c = 1/2 and the d = 1 transverse field Ising 
model estimate of x c — 1. Even without further anal- 
ysis these indicate the disordering effects of the strong 
frustration in the problem. 

Figure |^ also displays the excitation gap as computed 
via the modified Lanczos algorithm. Evident in the graph 
is the excellent agreement of the approximants and the 
exact diagonalization results over the region of consis- 
tency among the Pade approximants. 



1. Two Point Approximants and The Global Phase Portrait 



Although the region of validity of the analyzed gap 
series extends well into the domain of strong exchange 
coupling, the Pade approximants plotted in Figure |tell 
an ambiguous story. In particular, it is not clear from 
the displayed approximants whether or not the excita- 
tion gap vanishes for large J/Y and hence it is unclear 
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whether the triangular chain is ordered in the J 3> T 
limit. The reader should note though that even those 
Pade approximants that indicate a vanishing of the gap, 
do so for values of x in excess of 11. 

Biasing of the Pade approximants is one means of ex- 
tending the results of the gap series analysis to larger 
y values. A simple way to bias the-approximants is by 
means of "two-point" approximants. EZI In particular, one 
would hope that by biasing the value of the gap at y = 1 
the resulting approximants would be accurate over the 
entire range < y < 1. However, since it is precisely 
the value of the gap at y = 1 (let us denote that A(l)) 
that we know least well, such a procedure appears to beg 
the question. What we have done, then, is to construct 
two-point approximants (with various orders of numer- 
ators and denominator) for a range of A(l) values, and 
observed behavior. Two of the biased approximants, Pg 
and Pg , show a high degree of consistency with each other 
for a relatively small range of A(l) values. In fact, the 
range over which the maximum discrepancy between the 
two is less than 1% is confined between A(l) = .15 and 
A(l) = .2. For values of A(l) between .15 and .20, P 9 5 
and Pg 3 are in good enough agreement that the maxi- 
mum disagreement between the two over the entire y do- 
main falls below 1%. Outside the A(l) range indicated 
above, the maximum disagreement between Pg and Pg 
exceeds 1%. Figure |t] displays the biased approximants 
corresponding to the limits of A(l) discussed above. No 
matter what the choice of A(l), none of the other two 
point approximants are in agreement with each other or 
with Pg and Pg to an extent which approaches the level 
of agreement of Pg and Pg for A(l) values between .15 
and ..20. 

The biased Pg approximants displayed in the figure 
are in reasonable, though not perfect agreement with the 
Lanzcos gap values. The approximants shown are biased 
at A(l) = .15 and A(l) = .2, the bounds of the region 
in which P| and P| satisfy the 1% consistency criterion 
discussed above. The agreement is best for small values 
of y but deteriorates as y approaches 1. However, even 
at y = 1, averaging the two approximants yields a pre- 
diction for A(l) of .175, a value which differs from the 
Lanczos prediction of .132 (albeit without any attempt 
at scaling with system size) by only 30%. The approxi- 
mants echo the exact diagonalization results in indicating 
the absence of order for all values of J/T. 
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y = x/(x+1);x = J/T 

FIG. 7. Biased Pade approximant Pg plotted at both limits 
of 1% consistency in juxtaposition with Lanczos results for 
AT = 8. 



V. SUMMARY 

Results for the energy gap and cumulants from exact 
diagonalization of finite systems indicate that the trans- 
verse field, triangular chain Ising model is disordered, 
even though several analytic approaches (from both the 
small and large transverse-field limits) suggest that a 
particular ordered phase could exist for small transverse 
field. We note that the lowest lying excited state is con- 
sistent with the ordering analysis, indicating that the 
fluctuations generated by the ground state constraint at 
small T/J are too strong to allow order to set in. Hence 
we get "disorder by disorder" instead of "order by disor- 
der". Estimates of the gap obtained by direct Pade ap- 
proximants to 14th order perturbation expansions about 
the strong transverse field limit produce excellent results 
over a modest range of the Ising exchange coupling to 
transverse field ratio that exceeds the natural estimates 
for a critical value by a factor of 3 and do not indicate 
ordering for values of the ratio as big as 11. Globally rea- 
sonable results have been obtained by biasing the approx- 
imants and applying a consistency criterion and these 
indicate a lack of ordering at any value of J/T. This sug- 
gests that an application of this method to systems of 
greater experimental interest, such as the higher dimen- 
sional kagome and pyrochlore lattices, may be a viable 
approach to deducing the global phase behavior of such 
systems. We should note that in d > 1 we expect the 
series method to be increasingly competitive with exact 
diagonalization as the latter technique has to contend 
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with much stronger finite size effects at computationally 17 Carl M. Bender, Steven A. Orszag, Advanced Mathemati- 
feasible system sizes. caJ Methods for Scientists and Engineers, (McGraw - Hill, 

Inc., 1978). 
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